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ABSTRACT 

Possibihties are discussed for determining the top quark mass rrit from observations on the 
decay processes for top-antitop pairs produced in antiproton-proton colhsions at Tevatron 
energies, assuming that the t —>■ hW^ decay channel is dominant and much faster than 
hadronization. The final states ti — > bbfi^e^ provide the most striking signal, with little 
background, but they are rare (~ 2/81). If all three candidate events prove to be from tt, sua 
estimate follows for P{mt \ rate), the probability distribution for mt- The one reported con- 
figuration allows an independent estimate for P{mt \ ^'^e^ 2 jets). These two distributions 
are compatible and can be combined to give an mt estimate of about 122 GeV. Decay events 
"1 energetic lepton(/) + 4jets" should appear twelve times as often as "/x^e^ 2jets" events 
and can be analysed to give estimates for P{mt \ I4jets); this is illustrated for a fictitious 
event. There may be background from non-top events but suitable cuts on the data and our 
analysis procedure together reduce this to a low level. The rate observed for these events 
does not appear to be as large as this factor 12. Identification of either or both of the (bb) 
jets would be a great step forward, separating out the ''^^e'^W' and ''lbb2jets" events of 
importance with negligible background. We advocate an energetic approach to the analysis 
of individual events (whether 2,1 or no (6 or b) jets identified) on an event by event basis, 
with the hope of finding a subgroup of events with a common mass estimate.^] 



""^This paper was given by the first author (R.H.D.) on 19 July 1992 at the International School of 
Subnuclear Physics at ERICE. Its text has been updated to the end of 1992. 
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1 INTRODUCTION 



1.1 Existence of the top quark 

There is httle doubt today about the existence of a top quark t, the partner to the 
well-known and much-studied bottom quark b. The Standard Model, with SU(2) x U(l) 
for the Electroweak Interactions has had remarkable success, many of its parameters now 
being known to high accuracy. In particular, the weak-isospin component T^*^ for the b 
quark has now been determined empirically in a very direct way, from measurements of the 
forward-backward asymmetry in the process e+e" — >• bb and of the decay width T{Z^ — 66). 
Prom these data, Kane and Kolda [1] deduced the following T^'^ values for the b-quarks. 

T^\br.) = -0.49lH^ , T^'ibn) = 0.00l°;J° , (1.1) 

which are in accord with the Standard Model assignments of the b quark to the third quark- 
lepton family, for which the leptonic member t'^ is already known, as shown in Table 1. 

Table 1. The three families in the SU(2) X U(l) Standard Model 
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The fact that the predictions of the (spontaneously broken) SU(2) x U(l) symmetry have 
been so successful numerically provides a powerful argument that the third (quark-lepton) 
family must be complete. The absence of a T^^ = +i partner to the 6l quark would be 
such a gross violation of this symmetry that it would no longer be possible to understand the 
detailed and widespread agreement of the data with the predictions based on this electroweak 
symmetry. In short, there must exist a T^*^ ~ +| partner to the bottom quark b^, and it 
has been natural for this partner to be named the top quark t. 

Empirically, little else is known about the top quark, although its interactions are 
prescribed in form and magnitude from its place in the third quark-lepton family. However, 
we do know that its mass rrit is very much greater than the mass of the b quark, m;, ~ 5 GeV. 
A firm limit on rrit is provided by the fact that top-antitop pair creation is not observed in 
e~^e~ annihilation in the energy region of excitation; in other words, the threshold for tt 
production must be above Mz, giving the lower limit 

mt ^ Mz/2^i6GeV. 
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Indeed, the top mass is now believed to be above 91 GeV [2]. If top were lighter 
than this, tt pair production would have been observed already at the Tevatron, where 
proton-antiproton interactions are studied at cm. energy 1800 GeV, an energy far more 
than adequate for any reasonable value for m^. At the Tevatron, the limiting factor for top 
mass determination is the rate of events, rather than the available energy. We shall discuss 
this further in Sec. 2. 

The purpose of current top quark research is not just to demonstrate directly the 
existence of the top quark (which we believe already) nor to check that its interactions are in 
accord with SU(2) x U(l) symmetry, although the latter studies will provide a very fruitful 
field for research later on. The top quark will not lose its polarization by hadronization, as 
do the other heavy quarks, charm and bottom(c.f. Sec. 1.2), so that it will be possible to 
test all the detailed spin dependences of its interactions [3-5]. It may even turn out to be 
possible to use the top quark decay as a spin analyser for heavy particle and/or high energy 
processes which happen to give rise to top quarks. No, the paramount purpose of present 
top research is to determine the mass rrit of the top quark; how this measurement may be 
achieved will be discussed in Sec. 3 below. Only after this mass has been determined can we 
move on with the Standard Model, e.g. with the determination of the Higgs particle mass 
and with the testing of more extensive symmetry models which contain SU(2) x U(l) as 
a subgroup. Only then will we know what kind of accelerator we shall need for all these 
detailed studies. 

1.2 The top quark lifetime 

When rrit exceeds {M]y + m^) 85GeV, top quark physics becomes qualitatively 
different from the physics we have become accustomed to for the c and b quarks, as was first 
pointed out by Bigi [6]. The dominant decay mode for the top quark is then the two-body 
mode 

(a) t^W+ + b, (b) t^W' + b (1.2) 

As the mass rrit increases, this decay process becomes faster than hadronization. The top 
decay lifetime calculated for this mode [7] is plotted as function of rrit on Fig.l. For rrit 
appreciably lower than {M\,v + mj), the W-boson is virtual and its leptonic decay mode 
generates the overall mode t l^uih, whose partial lifetime falls like m^^ . In the transition 
region, as m* approaches and passes {Mw + rrib), the decay lifetime falls even more rapidly 
with rrit- For rrit well above this threshold, as holds for the physical situation, the decay 
lifetime falls more slowly, ultimately like m^^- The calculated partial width is given as 
function of rrit in Table 2. 



Table 2. Partial decay width for t bW^, as function of mt[29]. 



mt{GeV) 


100 120 140 160 180 200 


T{t bW+){MeV) 


88 298 612 1033 1572 2242 



[By way of contrast, we comment briefiy on the other heavy quarks Q = c and b. 
Their decay hfetimes are of order 10~^^ sec. In c and b jets, a polarized quark Q forms a 
meson {Qq) where q denotes a light quark q = {u,d,s), by picking up a light antiquark q 
from the vacuum. If this meson is in the ^5*0 ground state, all polarization carried by Q 
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is lost, since this state is spherically symmetric and cannot carry spin information. If the 
meson state (Qq) which is formed has non-zero spin, e.g. the "^Si state or orbitally-excited 
states Uke ^P2, it undergoes fast hadronic or electromagnetic transitions (since these conserve 
parity, the angular distributions do not depend on any initial polarization, although they can 
transfer polarization from the initial to the final state) until it reaches the ^5*0 state, which 
is the lowest state for the (Qq) system. Thus, all of the polarization information carried 
initially by the quark Q is lost, for the cases Q = (c, b)[8].] 

2 TOP MASS ESTIMATES FROM EMPIRICAL DATA 

2.1 Virtual top: radiative corrections 

It was first pointed out by Veltman [9] that, although the e+e^ energy at LEP is too 
low for the production of real top-antitop pairs, the existence of the top quark can still have 
a substantial effect on the predictions of the Standard Model through the "radiative correc- 
tions" generated by virtual top-quark loops and exchanges. Since the Standard Model is a 
renormalizable theory, these corrections are computable, at least in a perturbative approx- 
imation. Now that LEP experiments have measured with high accuracy many quantities 
related with the electroweak interactions, these measurements can be compared with the 
corrected theoretical predictions in order to draw some conclusions concerning the top quark 
and any other particles of high mass. Some of these measurements are the masses Mz and 
Mw, the total width Tz and some partial widths for decay, and the forward-backward 
asymmetries for e+e^ — > bb and e+e~ — > leptons in an energy range covering the peak, 
which arise from 7-Z° interference. In the minimal Standard Model, there are several other 
parameters also relevant, namely Mh, the Higgs particle mass, and q;5(M|), the QCD cou- 
pling strength evaluated at the Z^ mass. The latter can be deduced with fair accuracy from 
purely hadronic processes as well as from electroweak studies in the Z^ mass range. Given 
the recent LEP data and the theoretically computed expressions, it is then possible to lay out 
on an {rrit, Mh) plane, the regions consistent with these data. Quite a number of analyses 
have been carried out along these lines recently [10,11]. With the LEP data updated to July 
1992, EUis et al. [10] have given the value 

all data \mh free- f^t = 

124(27)GeV, (2.1) 

using as{Mz) = 0.118(8). Their analysis has some sensitivity to Mh', adopting by choice a 
rather low value for Mh, their analysis gives 

MH^Mz-.mt^ 132(25)GeV. (2.2) 

requiring a small increase for the optimum rrit. 

More elaborate models generally have more free parameters. The Minimal Super- 
symmetric Standard Model (MSSM) has aroused much interest recently, owing to its success 
in extrapolating the three coupling strengths, ai{n'^), at scale /x^ = (momentum transfer)^ 
and with i — (1,2,3) appropriate to the strong and electroweak sectors of the SU(3) x 
SU(2) X U(l) symmetry contained within a Minimal SU(5)SUSY, back to a common value 
at a GUT scale of order 10^^ GeV [12], a test that an earlier attempt based on SU(5) D 
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SU(3) X SU(2) X U(l) had failed [13]. Besides Mh and m^, this MSSM introduces five new 
parameters, itla and tan/3 for the Higgs sector, (^~, ""^o) the gluino mass and a common 
sfermion mass, respectively, and A^, a supersymmetry-breaking parameter. The discussion 
of the allowed regions for all of these parameters is naturally rather complicated. Ellis et al. 
[14] have discussed these dependences only for mt — 130 GeV, mainly to illustrate how an 
empirical value for would give much-needed information on the bounds which constrain 
the MSSM parameters. Most probably, given the number of new parameters introduced, the 
limits placed on mt in the MSSM will be much less restrictive than those for the minimal 
SM without supersymmetry. 

2.2 Total cross section for top-antitop production 

At the Tevatron energy of 900 GeV for proton and for antiproton, the (900+900) GeV 
pp interactions are rather like (300+300) GeV interactions between a valence quark q and 
a valence antiquark q. Creating a tt pair, each with (say) mt — 150 GeV , through the 
interaction process 

q + q^t + t (2.3) 



absorbs half of the initial energy, not leaving much energy for the creation of more particles. 
The energies of the residual (qq) and (qq) quarks from the proton and antiproton, respec- 
tively, lead through hadronization to particles which mostly go down the beam pipe. Thus, 
after the production process, 

p + p^t + t + X, (2.4) 



the t and t quarks decay according to (1.4 a,b), giving 

p + p^b + W- + b + W+ + X, (2.5) 



where X consists of (i) Xinv consisting of hadrons which go down the beam pipe, and (ii) 
Xyis, the other hadrons recorded by the detector. The bosons then decay, with lifetime 
3.1 X 10~^^s. Their simplest and most visible decay processes are 

(a) W+^l+ + ui, (b) W-^r + Pi, (2.6) 



for the three leptons I = e, jj, and r. 

The most striking events are those where the W'^ and W~ decays lead to two charged 
leptons from different families without any X^^s, a typical final state being 

b + b + ii^ + e^ + (ui, + Pe + Xinv), (2.7) 

where the b and b quarks hadronize to give corresponding jets j{b) and j{b). There may also 
be some secondary jets, emitted from the initial quarks or by the quarks heading for the 
beam-pipe, or arising from the development of the b and b jets, but these will generally have 
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relatively low energy, so that these {^^e'^bb) final states may be expected to have rather 
simple structure. Another attractive feature of these (fi'^e^bb) events is that they may be 
expected to be rather free of background, since n and e belong to different families. There is 
no particle Y known whose decay leads to the p,e or //e configuration without corresponding 
neutrinos ^'e(or Pg) a-nd ^/^(or i^^). Indeed, lepton conservation is believed to hold separately 
for muonic and electronic lepton numbers: strong limits have been placed on any violation of 
this conservation law, through searches for reactions of the type — >■ e~p and for decays 
of the type Y — > /i^e^X, especially for Y^ = and Y^ = . 

The simplest pp process leading to final states with {jJi^e^) is VT-pair production 

p + p^W- + W+ + n jets, (2.8) 

followed by the decays (2.6). This process involves two weak vertices whereas the process 
(2.4) has none; provided that the energy available is reasonably large compared with 2{Myy + 
nib), the process (2.8) with n = 2 may be expected to have a rate smaller by two factors 
of (aem/'^sin'^Ow) than that for (2.5), if the n — 2 jets were 6-jets. However, the process 
(2.8) will be in fact dominated by events with light quark jets in large numbers, and it is 
difficult to estimate how much background (2.8) will produce for the tt events (2.4) leading 
to final state (2.5). It is a question of what fraction of these non-tt events (2.8) having 
n = 2 light quarks and surviving the transverse momentum cuts (c/. Sees. 3 and 4) can be 
adequately fitted to the interpretation (2.7). The only further remark we can make is that, 
if a microvertex detector is used, so that it is possible to add the requirement that two of 
the jets in (2.8) be identified as b jets, the background rate from (2.8) would be of order 
{Q;em/(2as<sm^^^vv')}^ ~ 10^^ times the signal for (2.7), following the top-pair production 
(2.4). In our further discussion here and in Sec. 3, we shall neglect background to the 
process (2.5) arising from the direct production of W pairs. However, the validity of this 
neglect is under discussion at present. 

The cross section for tt production in pp annihilations at the Tevatron energy of 
900 GeV through reaction (2.4) has been calculated by Eichten [15], as function of m^. At 
Tevatron energies and for such large rrit values ( ~ 120 GeV) as we are led to consider 
here, it production is dominantly due to the quark-antiquark process (2.3); production by 
gluon-gluon processes is small in comparison. 

The bosons have hadronic modes, in addition to the leptonic modes (2.6). These 
are due to the W^-coupling with the quark weak-isospin charged current. Since we are inter- 
ested here only in total rate, we limit our discussion to the couplings 

W-^ ^ (a) u + d, (b) c + s, (2.9) 

given by diagonal terms of the Cabibbo-Kobayashi-Maskawa matrix. Each of these couphngs 
leads to two final quark jets. The third coupling is W~^{t + b) whose threshold lies far 
above M^, so that it contributes only to virtual processes, such as the radiative corrections 
mentioned above. 

Recognizing colour, there are 6 couplings in all, each with the same strength as 
each W'^ — > I'^ui. Empirically, the branching ratio for each (lui) final state is 10.5(9)%, 
which agrees with the naive expectation of 1/9, from the channels just enumerated. If both 
and W~ in (2.5) decay hadronically, the final state will have 6 energetic quark jets, a 
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configuration which it may be difficuh to disentangle. If only one W decays hadronically, the 
final state will generally have one energetic lepton and four quark jets. The analysis of this 
final state is generally possible and we shall discuss this in Sec 3. The dilepton events have 
been described in part just above. The decays (2.6) include the possibility of emission; 
this is far more difficult to deal with experimentally than is or n^. For this reason, we 
shall not refer explicitly to the possibility of emission again. On the other hand, final 
electrons or muons can be easily recognized and distinguished. Among the dilepton final 
states, the (e+e^) and cases need special attention and much caution, since there 

are many other mechanisms which can give rise to electron or muon pairs, one example being 
the production and decay of a Z° boson, and which will still require special cuts for their 
exclusion. In fact, no candidate event of this type has yet been reported, although such 
events (e"'"e~ and must have the same rate as the (/^^e^) events. 

2.3 Top mass estimate from production rate 

In this section, we confine attention to the {ii^e^) events. The expected number iV^g of 
these events is plotted as function of rrit in Fig. 2 for an integrated luminosity IL = 30p6~^, 
based on the estimates by Crane [16] who used Eichten's calculation of the cross section 
for top-antitop pair production [15] and took due account of the efficiency k for the CDF 
detector as a function of the event location and configuration. For mass m^, we have 

N^, = K-IL-a^e{nit). (2.10) 

where cr^g denotes the cross section for (n'^e^) events in pp coUisions at 1800 GeV cm. en- 
ergy. The probability of producing n such events, given m^, is 

P{n I mt) = ^^^^^exp{-N,e{mt)) (2.11) 

77/. 

Given IL, (J^e("^t) a-nd the observation of n events, the Bayesian probability distribution for 
rrit is then 

M (m V" r N (ni)"' 

P{mt I n) = ^^Y^^ exp{-N^,{mt))/{ / dm ^ \ ' exp{-N^,{m)}, (2.12) 

77/. J 77/. 

where the integral is taken over all possible values m allowed for the top quark t by all the 
conservation laws. This probability may be reduced to the form 

P{mt I n) = c{n){K ■ IL ■ (T^e(mt))".exp(-«; • IL ■ a^eim)), (2.13) 

c{n) being an n-dependent normalization factor. This distribution peaks at the mt value for 
which 

(T^,eimt) = nlK- IL . (2.14) 

One good (/i~e^) candidate event has already been published by the CDF collaboration 
[2] and has been discussed [17] in some detail (c/. Sec. 3 below). A second (/xe) candidate 
was shown by the CDF collaboration in their report given at the November 1992 Chicago 
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Meeting of the Division of Particles and Fields of the American Physical Society, although 
no measurement details were released. It was well known at that meeting that the DO 
collaboration also had their first (/le) candidate. | Although the integrated luminosities IL 
are not known to us precisely, a value of about 20 pb~^ for CDF (including IL=4.7 pb"^ 
from their 1989 paper) and 10 pb^^ for DO would appear plausible estimates, at least of the 
right order of magnitude. The mean detector efficiency can be deduced from a comparison 
of Crane's rates with Eichten's total tt production cross sections. 

On the assumption that these three (/ie) candidates do stem from top-antitop 
production, and that the integrated luminosity up to November 1992 was about 30 pb"^, the 
probability distribution for nit is shown on Fig. 3. Its peak is at 120 GeV, the one-deviation 
limits being 109 and 135 GeV. Since the curve for iV^e(^t) shown in Fig.2 is a steeply falling 
function of rrit, the peak value thus determined for rrit is not strongly dependent on our 
estimate for IL, nor on the number of /le events. For {IL,n) = (40p6~^,3), the peak value 
is mt = 127.5; for (30p6~\4), it is at 114 GeV, and for {30pb-\2), it is at 129 GeV. 



3 ANALYSIS OF DILEPTQN EVENTS 

3.1 Kinematics of top decay sequence t — > bW~^ — > bJ^iyi 

Consider first the kinematics of top decay to bW~^, followed by decay to lui, in 
any frame. Energy-momentum conservation gives 

t = b + l + i^i. (3.1) 



where a denotes the energy-momentum four-vector of the particle a. Since / and ui are decay 
products of W, we have 

M^ = {l + uiY = {t-b)\ (3.2) 

and, since the neutrino has zero mass, 
= (^t-b-lf_ _ (3.3a) 

= (t-bf -2l.{t + b) + l\ (3.36) 
where a.b denotes the scalar product of the four-vectors a and b and = a.a. For / = e or 
/i, we can neglect the lepton mass, i.e. l'^ = mf = 0. Using (3.2) in (3.3b), we obtain the result 

l.t = Lb + (3.4) 

Evaluating the r.h.s. in the lab. frame and l.h.s. in the top rest-frame, we deduce that 
Eu = ilb + M^/2)/mt, (3.5) 

giving Eit, the lepton energy in the top rest-frame, in terms of lab. measurements for I and 
b. 



^This event was shown at a Fermilab seminar in February 1993, but without full details. 
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It is useful to derive here an inequality for mj. From (3.1) and (3.2), we have 



mf = {b + l + uiy _ (3.6a) 

= ml + + 2b.l + 2h.vi (3.6fe) 

Evaluate the invariant products h.l and h.vi in the W rest frame, where Ei = I and Ey = v 

are each My/ and the momenta / and vi_ are opposite. Denoting by 9 the angle between I 
and h in this frame, their product takes the form 

hlh.vi = {Eh-bcose){Eb + bcose)(Mw/2y (3.7a) 

= {ml + b^sin^9){M^/4:) (3.76) 

> mlM^/A. (3.7c) 
Using this inequality for h.ui in Eq.(3.6b), leads to the inequality 

m^^ > ml + M^_+2b.l + mlM^/2b.l (3.8a) 

- {ml + 26./) (M^ + 26./) /26./ (3.86) 
as given in the Appendix of ref.[5]. 

3.2 A geometric construction 

Now consider the kinematics in the lab. frame, starting with the 3-momenta b and /. 
The top momentum t is constrained by two relations, 

{t-bf = {E-E,f-M^, (3.9) 



from eq.(3.2) and 

{t-b-lf = {E-E,-Eif, 



(3.10) 



from eq.(3.3a), where E has not been constrained to the value Et because the value of rrit 
is still left open. For given E, eqs. (3.9) and (3.10) constrain t to lie on the intersection of 
two spheres, one centred on the point 6 and the other centred on (6 + /). Their intersection 
is a circle on a plane perpendicular to the line BL joining these two centres and centred on 
this line. As E varies, the centre of this circle moves along the line BL at a rate linear in E, 
and the square of its radius also increases linearly with E; in short, as E varies, this circle 
traces out a paraboloid with axis BL, as described in refs. [8,17]. All points t satisfying the 
constraints (3.2) and (3.3a) for given 6 and /, for all possible rrit values, lie on this paraboloid. 
The points of interest to us are those for a definite nit, still to be determined. These points 
lie on an ellipse formed by the intersection of the paraboloid by a plane whose normal lies 
in the plane OBL and makes an angle a with the axis BL, where 

tana ^bi/iEb-bs), (3.11) 

is independent of nit, and 6i and 63 are defined on Fig.4. For a specified rrit, the possible 
vectors t are given by OT, as T moves around this ellipse, and the vector LT gives the 
corresponding neutrino momentum z^. It is apparent that the top quark energy Et = \/(t7T.| + 
t^) must lie between two limits Emin{'^t) ^ Et{mt) < Emaxi^t)- 

As rrit increases, the plane PTQ moves upwards, with a constant normal since the angle 
(7 (given by (3.11)) does not depend on rrit- The ellipse retains the same eccentricity but 
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increases in linear dimension. As rrit decreases, the plane moves downward and the ellipse 
shrinks until for a limiting value nit = '>ti* the plane becomes tangent to the paraboloid and 
the ellipse is reduced to a point. For rrit < m*, there is no solution for t. This value 
is the lowest limit for an rrit consistent with the momenta (6,1); m* is, of course, the value 
given by the square root of the expression (3.8b). The projection of the ellipse for rrit on 
a plane perpendicular to the vector Z is a circle, and the centre of the ellipse projects onto 
the centre of the circle. As rrit varies, the centre of the ellipse moves parallel to the vector 
BL = I and is distant from it by amount {M^/2Ei)tana in the plane of b and /. 

It is worth noting here that any configuration {b,T) can be fitted by a sufficiently 
large rrit- However, a very large rrit means that LT is very large, but this is the neutrino 
momentum z^, of which we have no direct knowledge. It means also that OT must be 
correspondingly large, and this is the top momentum t. These two very large momenta 
would be compensating each other, differing by the observed momentum (6 + /), and this is 
not generally a plausible interpretation of the top decay event observed. 

For top-antitop production, with lepton emission in both decays, two paraboloids are 
to be constructed, one from (6, /) for t and one from (6, /) for t, both in the same laboratory 
frame. One dilepton event has already been described in detail in the literature [17] from the 
CDF experiment in 1988 as a possible t production event. We shall refer to it as CDF- 1. Its 
measurements are given here in Table 3, the z-axis being along the beam. There is a third 
lepton in CDF-1, a //"^ with quite high energy but with rather small transverse energy. It 
allows interpretation as a secondary muon, emitted from B meson decay in the 6-jet following 
the i hW'^ decay of the hypothetical antitop quark produced in this pp interaction. It has 
the right charge sign and a low transverse energy compatible with the sequence h —>■ cW^, 
W'^ — > i^^Pij,, and it travels in the same general direction as an energetic hadronic jet. These 
features strongly suggest that the and the hadronic jet are to be taken together, as the 
components of a 6 jet. We expect that the other jet is a b quark jet but there is no evidence 
to demonstrate this. 

3.3 Top mass estimation from dileptonic decays i — > fl^e^ + 2 jets 

If we now assume that CDF-1 represents top-antitop production, with the assignments 
given in Table 1, the two paraboloids then constructed are those shown in Fig. 5, projected 
onto the (y,z) plane. The eUipses representing t decay are clearly seen; the i-eUipses happen 
to be seen edge-on in this projection, since the momentum vectors b_ and I almost lie in this 
plane, their components and being small. These two sets of ellipses are shown for rut 
ranging from 115 to 305 GeV, in 10 GeV intervals. The quark- ant iquark interaction which 
gave rise to this event was strongly asymmetric, since the two paraboloids lie almost entirely 
within the same hemisphere. 

It had been assumed generally, we think (but see ref.[19]), that it would not be possible 
to analyse such an event in terms of top-antitop pair creation since the final state would then 
include at least two energetic neutrinos, z/g and z/^. [CDF-1 actually has a third neutrino, of 
type and with unknown energy, from the decay mentioned above. It is assumed here that 
this energy is small although its associated has quite high energy. This assumption 
needs checking for internal consistency]. Of course, we do not know whether this t assumption 
for event CDF-1 is correct, but if it is, then it is possible to deduce something probabilistic 
about rrit- 

Table 3. Input data for two t candidates, one real (CDF-1) and the other fictitious (F^e), 
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both of the form [ji^e^hh). 





Event CDF-1 


Event F,,^ 


GeV/c 


e+ jet(6) e- jet(6)* 


e" jet(6) //+ jet(6) 


Px 

Py 

Pz 


-21.2 18.7 0.6 -0.4 

23.6 -6.25 -43.7 16.8 
-28.55 25.25 16.8 -98.9 
*includes a low energy 


22.5 66.3 8.0 -47.7 

-5.45 86.4 20.0 9.5 
-21.5 -12.0 94.0 27.5 



There is one further input. In the simplest parton picture of the pp collision, [20], each 
of the partons q and q has no transverse component of momentum. It follows that the net 
transverse momentum of the final t and t quarks must be zero, so that we have the two 
further constraints: 

ty + t^ = 0. (3.12) 

where the suffix T denotes the component of the vector transverse to the incident proton 
beam. What is the consequence of this condition? We can see this most simply by ro- 
tating the paraboloid by 180° about the z-axis. If we view the two resulting paraboloids 
along the z-axis, i.e. project them on to the XY plane, the constraint equations (3.12) for 
an assumed rut correspond to a crossing of the two ellipses for this mt, since (3.12) then reads 

-tT^tT- (3.13) 



These crossings are shown for a number of rut values on Fig.6, for the event CDF-1. For 
rrit < 110.2 GeV, the ellipses do not cross (for rrit < 100 GeV, they do not even exist) nor 
do they cross for m,t > 410 GeV (which is an unreasonably large rrit, as is explained below). 
For rrit between these limits, there are two or four crossing points, each giving a solution for 
t and t (and hence for i/g and separately) . The next step is to assign a probability to each 
configuration. 

From the z-componcnts of a solution for the momenta (t,t), the momenta xP and xP 
of the parton and antiparton, where P denotes the total momentum of p and p in the Lab. 
frame, can be deduced, giving 

x^{Et- E-t + {tL + fL))/2P, (3.14a) 



and 

x^{Et- E-t - {tL + iL))/2P, (3.146) 

where the suffix L denotes the component of the vector along the incident proton beam. 
The structure functions for the proton and the antiproton are the same and are well-known 
[21]. For the proton there are three of these, F{x) for the quark partons, F{x) for the 
antiquark partons and Fg{x) for the gluons, and they depend on (momentum transfer)^, for 
which we have taken the value (m^)^. For the antiproton, the functions F{x)andF{x) are 
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to be interchanged and x replaced by x. For the values of {t,t) relevant to the Tevatron 
experiments, the relevant values x and x are quite large, being typically 0.1 or larger. For 
such values of x and x, F is small and Fg is even smaller. We should add here a note that, 
some way below the upper limit allowed for by the kinematics, the values of x and x 
required by eqs.(3.14) may exceed unity, as rrit increases; of course, the expressions (3.14) 
are then not physically meaningful solutions. Thus, the condition that x and x must be less 
than unity imposes an upper limit on the magnitude of nit. 

Thus the it production rate at the Tevatron is dominated by the quark- ant iquark 
collisions, which are proportional to F{x)F{x). We have used for F set 2 of the Duke and 
Owens [22] structure functions for (u+d), neglecting F and Fg. Even F(x) falls to zero quite 
rapidly, like (1 — x)^, as x increases above 0.5. The rate is also proportional to the differ- 
ential cross section for parton — antiparton tt, summed over all of the parton-antiparton 
collisions. Including for completeness the i = {1,2,3} = {{qq), (qq), {gg)} initial states, we 
find for the rate factor, 

R{x,x) ^j:F,ix)F,ix)^^%^, (3.15) 

where Fi{x) now stands for F{x), F2{x) for F{x) and F3{x) for Fg{x), and s — xxs, and 
t — mf — {Ft — tijx^s are the differential cross section variables. 

Finally, we have to consider the spectrum of the lepton in the top rest-frame (and of the 
antilepton in the antitop rest frame). With the Standard Model couplings W ■ {^^{1 + 75)6) 
and W ■ {I'yil + 75)t'0> there is a large forward-backward asymmetry in the W — > lui decay 
relative to the b-momentum in the W rest-frame. Since there is a one-to-one relationship 
between this angle and the energy Eff we can express this angular distribution as an energy 
distribution: 

-— = GlM^Ejtimf -ml- 2mtFit) / {A^nitTw). (3.16) 

If we neglect {nib/ nit) << 1) and normalize the distribution to unity, we have the probability 
distribution 

P{Ert) = {24/m't)Eit{l - 2Eit/mt)dEit (3.17) 

where Fjt runs from to mt/2 in this approximation. 

Taking all these factors together, the total probability for reaching the observed con- 
figuration is proportional to 

Pill I mt) ■ ^^^^ ■ P{Eu I mt)PiEu \ m,), (3.18) 

j i "'''i 

where the sum i is over (qq), (qq), (gg), and the sum j is over all of the ellipse-crossing points 
for the mt value under consideration. As discussed here and in ref.[18], P{t,t \ mt) is a dis- 
continuous function of mt, the ellipses being lines of zero width which either cross or don't 
cross. In reality, these lines have finite width, partly because top quarks and W-bosons have 
finite widths, but more importantly because of statistical uncertainties in jet development 
and in measurements of the particles in the event, and because of a distribution of transverse 
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momenta for the partons and the antipartons in the incident proton and antiproton. In order 
to estimate rrit from the data on CDF-1, assuming it to be due to top- antitop production 
and decay, we appeal to Bayes Theorem, which gives the probabihty distribution as: 

P{'>TT't \t,i)— P{t, 1 1 Mmt)/{ J dmP{t, i I m)$(m)}, (3.19) 

where $(m) is the a priori probabihty that the top mass is m and the integral is taken 
over all possible values of m. With $(m) = 1, i.e. for the CDF-1 event shape alone, the 
probability P{mt \ t,i) is plotted in Fig. 7(a), as taken from ref.[18]. Its peak lies at about 
130 GeV, and it is broad and flat. 

3.4 More realistic top mass estimation. 

Although the above analysis is geometrically simple and easy to comprehend, we have 
come to reahse, from the analysis of the "l-|-4jets" events in Sec. 4 below, that it is really 

necessary to allow (i) for the transverse momenta of the initial parton-antiparton system, 
and also (ii) for the uncertainties of the b and b jet energies. In this Section, for simplicity, 
we shall consider only the former of these; the inclusion of the latter is discussed below, in 
Sec.4.2. 

Consider first the t-ellipse for given b and I, for assumed top mass m. Divide its 
boundary by N points by dividing the ellipse of Fig.4 into equal azimuthal slices with respect 
to the axis CeN, where Ce denotes the centre of the ellipse and CeN intersects BL at angle 
a given by expression (3.11) and is normal to the plane of the ellipse. At each one of these 
N points (a set labelled by 7) there is a definite momentum ^(7). Next, consider the t-ellipse 
for b and I, for the same value of m, following the same procedure for the N points (this set 
labelled by 7), each point leading to a definite momentum ^(7). We now soften the condition 
(3.12) by the introduction of a weighting factor Tp of finite range, thus obtaining in place of 
(3.18): 

P{t,i\ m) = Y.{P{t^ I m)Tp{t^T + hT)P{h I 
77 

where T'p(t) is the function 

which is a representation of the two-dimensional 5-function 



(3.20) 



(3.21) 



(3.22) 



as p — 0. In practice, we will use the sum (3.20) over the set of points (7,7), each having 
180 points, so that the slices are of angle 2°. The contribution of each point to the sum is 
evaluated and summed. The largest individual contributions are naturally those from the 
vicinity of the crossing points of Fig. 6. The results of this calculation are shown in Fig. 7(a) 
for p = O.lmt and 0.025mt.We see there that the forms obtained differ considerably from 
that obtained with the condition (3.12), and further, that there is no sign that the earlier 
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result is recovered in the limit p — 0. However, it is easy to understand the change in form. 
Consider the case of m = 125 GeV. We no longer have single large contributions coming 
only from the crossing points C and D, but find a large number of smaller contributions all 
the way from C to D; this situation persists as the mass falls below 125 GeV, even after 
there arc no crossings corresponding to C and D. Since the ellipses become smaller as m 
falls, the finite range gains increasing relative importance and it is not surprising that the 
non-crossing contributions become dominant there, giving greater weight to low m values 
and thus distorting the m-distribution away from the form which we presented in ref.[18]. 

The dashed line on Fig.7(a) shows the m-probabihty of the single point having the 
greatest probability for the m considered. This point naturally lies close to the cross-over 
point for that m, if there is one, so that this entry must be closely related with the value 
obtained with the procedure based on (3.12). This curve, based on these entries of maxima, 
does have the slow fall-off for increasingly large m which we noted in ref.[18] and show here 
in Fig. 7(a). The change introduced by the use of (3.20) arises from the fact that there are a 
large number of points (7,7) which cannot contribute to the sum P{t,t \ m) when (3.12) is 
required but which dominate its sum for low m when (3.20) is used, so depressing the value 
of the most likely value for m. The probability distribution we now adopt is that shown in 
Fig. 7(a) for p = O.lrrit. 

3.5 Compatibility and combining of top mass probability distributions. 

The probability functions P(mt | iV^e) and P{mt \ t,t) are independent and at present 
compatible. They can be combined to give an overall probability function from the data 
publically available today: 

P{m I data) = P{m \ N^e = 3) • P(m | t,t:CDF - 1), (3.23) 

shown on Fig. 7(b). Its peak is at 121 GeV, the one-standard-deviation limits being 114.5 and 
130 GeV. Of course, these remarks are very tentative since they are based on the assumption 
that the three (//^e^) events are due to top production and decay; all three might be due to 
background. 

At present, we do not know whether or not rrit = 121(+9, — 6.5)G'ey is the top 
mass. However, the hope for the future is that a substantial number of (/x^e^) events will be 
observed at the Tevatron. An rrit distribution can be determined for each of them, as done 
for CDF-1 here. These distributions may turn out to vary randomly from event to event; if 
so, then these events could not have anything to do with top-antitop production and decay. 
More likely, a large fraction of them (and perhaps all) may peak at a definite rrit value and 
we shall be able to conclude that they form a well-defined group of events due most probably 
to ti production and decay. The total cross section for this group of events will provide a 
quantitative test on this interpretation for them. 

In order to illustrate one difficulty in following this line of thought, wc shall close this 
Section by discussing a fictitious event which we shall refer to as F^e- Its input momenta 
are given in Table 4 and it has been subjected to the same analysis, with p = O.lm^, as was 
applied to the real event CDF-1. Evaluating the limit (3.8b), the antitop momenta lead to 
{mt)min — 95.36*6^, while the top momenta lead to {mt)min — iM.BGeV. Since the latter 
is the stronger limit, we conclude that, if F^^ were data from a real tt event, then the top 
mass could not lie lower than 134.5 GeV. The nit distribution for F^e is shown on Fig.8. 
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It peaks at 154 GeV, with one- standard-error limits at 144 and 168 GcV and it is almost 
incompatible with the rrit distribution given in Fig. 7(a) for the real event CDF-1. We would 
have to conclude that either one (or both) of these events is not the result of top-antitop 
production and decay, or that there has been a large statistical fluctuation. In either case, 
the situation would be unsatisfactory and we would have to await further /xe events. At the 
moment, as mentioned in Sec. 2(b), it is known that there do exist two further events of this 
kind and we all look forward to the release of the detailed data from them. 



4 ANALYSIS OF "LEPTQN AND FOUR JETS" EVENTS 

4.1 Expected rate and nature of the events. 

As mentioned following eq. (2.9), the hadronic decay modes — > udandcs have a 
net rate about six times greater than either of the decay modes W'^ — > e+^e or /x"*"!/^. From 
this remark it follows that tt decay leading to the flnal charged particles 

e"^ or iJ,^)b) + b{udorcs), (4.1a) 

and 

{l~ {= e~ or fj,~)b) + b{ud or cs) , (4.16) 

are more frequent than the {/i^e^) dilepton modes 

(/x+6) + (e-6), (4.2a) 



and 

(/.-6) + (e+6), (4.26) 

by a factor of twelve. Here, we have confined attention to the diagonal elements of the 
Cabibbo-Kobayashi-Maskawa matrix and then only for A = (Wolfenstein's parametrization 
[28]), since this is a satisfactory approximation for the present. If the three (//^e^) events 
mentioned in Sec. 2(b) were all due to top-antitop pairs, then we would expect about another 
36 pairs, each giving one energetic lepton with high accompanied by one b~jet, one b—jet 
and two other jets, as specified by the possibilities given in (4.1). Up to the present, it has 
not been possible to identify which are the 6 and 6 jets, except by the observation of a 
low-energy muon resulting from the secondary decay b ^ l~Pic or b ^ I'^vi, as was the case 
in the CDF-1 event, where there was a low-energy ^+ (see Table 3). The secondary lepton 
from 6 or 6 decay has branching ratio about 21%, so that one should appear in ^ 42% of tt 
events; the appearance of two secondary leptons has a rate an order of magnitude lower. In 
the present run, CDF is equipped with a microvertex detector, which will greatly increase 
its efficiency for identification of the 6 and 6 jets, because of the relatively long lifetime of 
the b quarks (pa 10~^^sec), and so restrict the fitting of the final state to the tt hypothesis 
in a very helpful way. 
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Table 4. Input momenta for a fictitious tt event (F;) of the form {I bbq'q), where (g', q) — 

(s, c) or {d, u). 



GeV/c 


e 


jet(e) 


jet(V') 


jet(0) 


jet(r/) 


Px 


11.0 


42.0 


0.7 


0.0 


7.0 


Py 


-19.0 


-23.5 


23.5 


-25.5 


16.0 


Pz 


-7.0 


25.0 


-56.0 


23.0 


72.0 



4.2 Their analysis and top mass probability distributions. 

Here we discuss the analysis of the (/ -|- 4jei) events, when there is no information as 
to which are the h and h jets. What we say below will assume the lepton to be the l'^] if it 
is then the terms top and antitop are to be interchanged in what follows. 

We start by labelling the jets according to their momenta, using the symbols (0, V', e and rj), 
as illustrated in Fig. 9 and in Table 4. The mass Mw — 80.2{3)GeV and the width — 
2. 1(1) Gey are now rather well known for the W bosons; we shall therefore impose this value 
of Mw on the event, while neglecting the width for convenience, since its effects are of 
secondary importance. Two of the jets are selected, say e and rj, and identified as 1 and 2, 
with energies Ei and E2 and separation angle 6*12. We then have 

M^^ml + ml + 2{EiE2 - pmcosOu) ^ 2EiE2{l - cosOu). (4.3) 

Each jet energy Ej has a probability distribution Qi{Ei), partly because jet development is a 

stochastic process and partly because of uncertainties in its determination from the observa- 
tions. Eq.(4.3) fixes the product E1E2, so that the allowed values correspond to a hyperbola 
on the {El, E2) plane. A''i2 points are chosen on that hyperbola and the probability assigned 
to each of these points is deduced from the integral 

g,,(l, 2) = y dEidE2Qe{Ei)Qr,iE2)5{EiE2 - ^^73^^)- (4-4) 

A third jet, say 0, is selected and labelled 3; A^3 points are chosen over the energy range 
(£■3 ± xds), where cts is the assigned uncertainty at the one-standard-error level and x is 
chosen suitably (usually x = 2). Thus, we have chosen A^3 ■ A''i2 energy values (a set labelled 
a); for each member of this set, an antitop momentum t{a) and mass m{a) are deduced, 
with an assigned probability Q{a), which is the product of Q12, computed from (4.4), and 
Q4>{E3). 

We now turn to / and the fourth jet labelled 4. This jet is tacitly (and necessarily) 
assumed to be a b jet, although this identification is not used here. We take N4 energy val- 
ues (a set labelled /3) over the range £^4 ± X(J4, each with probability Qij}{P). For each point 
{a, (3), we determine the ellipse for mass ■m{a), based on the vectors I and b^, and divide 
its boundary by A^; points (a set labelled 7), as we did in Sec.3.4 by taking equal azimuthal 
slices of the ellipse of Fig. 4 with respect to the axis CeN, where Ce denotes the centre of the 
ellipse and C^N intersects BL at the angle a given by expression (3.11) and is normal to the 
plane of the ellipse. At each one of these A^; points, there is a definite momentum t{a, (3, 7) 
for the top quark. Next, we ask for a match between t and t, i.e. that 
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t_{a,p,j)T + t{a)T = 



(4.5) 



if the incident parton-antiparton system is required to have no transverse momentum. How- 
ever, in general, this equation will not be satisfied, except at isolated points. As above, in 
Sec. 3.4, we give up the equality (4.5), allow initial transverse momentum and introduce 
the probability distribution (3.21). where here r is given by the l.h.s. of cq.(4.5), adopt- 
ing as before the value p = O.lm^. Just as for the dilepton modes in Sec. 3, a lepton factor 
(3.17) is necessary, which we denote here by P^{Eif, 4). For the accompanying decay sequence 
t — > bW~, W~ — > cs or ud, the light quark s or d will have the same strong backward/forward 
asymmetry as holds for the lepton in W iPi, but since we have (in general) no means to 
distinguish the d jet from the u jet (or the s jet from the c jet), we have to average over the 
two final jets and this removes the backward/forward asymmetry. For simplicity, we have 
taken the decays i — > bcs and bud to be isotropic. For the production of i and t quarks 
with momenta t{a) and t(Q;,/5, 7), through the process qq — > it, there is also a rate factor 
R{x{a, P,j),x{a)) necessary, given by (3.15) with x and x defined by eqs.(3.14). Thus, for 
the labelling (3,4,1,2) of the four jets (0, ip, e, rj) and for each of the number (iV/ ■ iVi2 ■ ■ N4) 
of points defined above, we have calculated a mass m(a), the two momenta t{a) and t{a, [3, 7) 
and an associated probability for reaching this event configuration, given by the product 

P{event\m{a); 3, 4, 1, 2) = Q,^{1, 2)Q^{'i)Q^{A)T{a, P, i)P^{Ej,, A)R{x{a, P, j),x{a)). (4.7) 

Finally, for this particular assignment of the labels (1,2,3,4), the net probability for this 
event configuration when ni{a) lies in the interval {m,m + Sm) is obtained by summing 
these probabilities (4.7) for all of these {Ni ■ N12 ■ ■ points for which their m{a) lies in 
this interval. 

When none of the four jets is identified, it is necessary to consider all possible iden- 
tifications for them. The final expression for the net probability that the observed event 
configuration could occur for top quark mass m is given by the sum of (4.7) over all permu- 
tations of the labels (1,2,3,4): 

P{event \m) = J2 P{event \ m; 1,2,3, A). (4.8) 

Perm. 

As in Sec. 3 above, we now use Bayes' Theorem to obtain the desired expression for the 
probability that the top quark mass is m, given the data on this event and an a priori prob- 
ability $(m) from other information (such as that given by the tt production rate observed; 
cf.Sec.2.2 above), as follows 

P{m I this event + earlier info.) — c ■ P{event \ m)$(m), (4-9) 

where c is a normalization constant, independent of m. 

Of course, most of the identifications in the sum (4.8) are necessarily inappropriate 
but we have no means, in general, to know which of them correspond to the correct physical 
interpretation. It is our optimistic expectation that incorrect identifications will generally 
lead to negligible contributions to the net probability for top mass m, although it can cer- 
tainly happen that several different interpretations have comparable probabilities. The merit 
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of this method of analysing the data on an event is that it is systematic, so that no possible 
interpretations can be overlooked, and that it is quantitative, assigning a numerical value 
for the relative probability for different interpretations. It is not an elegant procedure but it 
does not involve multivariate searches for the maximum likelihood. 

We may illustrate the procedure by discussing a fictitious event of the type (l + Ajets), 
with final state momenta as given in Table 4, which we shall refer to as Fi. The probability 
distribution P{m) calculated for Fi is shown on Fig. 8. This distribution peaks at 135 GeV, 
with one-standard-deviation limits at ±6GeV; it has quite a strong overlap with those for 
CDF-1 and it overlaps with that for F^e no more than that for CDF-1 does. 

4.3 Simulated data and its analysis. 

To test the method more extensively, we have considered two sets of simulated data: 

(a) A Toy Model. 

This is described in ref.[23]. Adopting a mass of 140 GeV, the vectors t and t were 
chosen in a random way, and their decay configurations follow the leptonic sequence for t 
discussed in Sec. 3 and an isotropic 3 jet sequence for t. These events were processed by the 
procedure described above, including allowance for the Gaussian probabilities appropriate 
to the CDF determination of jet energies from the data on each jet and for all permutations 
of the assignments of (1,2,3,4) to the jets. In this analysis, the same cuts were made on the 
simulated data as are routinely applied by CDF to their real data. These cuts are as follows: 

(i) pt > 15GeV for each jet, 

(ii) EiT > 20GeV for the lepton, 

(iii) missing transverse > 35GeV, 

(iv) pseudorapidity lower than 2.44 in magnitude for all four jets. 

The P{m \ event) distributions for three "toy events" chosen at random from 1000 
simulated events are shown in Fig. 10. In two of them. Figs. 10(a) and 10(b), the peak of 
the distribution hes close to 140 GeV, the input value although with much spreading due 
to wrong jet combinations. However the third, Fig. 10(c), shows that individual events can 
deviate widely from our simple expectations, for its peak is at 150 GeV while it has a marked 
dip in the vicinity of 140 GeV. In Fig. 11, mean probability distributions for 1000 "toy model" 
events generated for rrit = 140 GeV, are shown for two cases. For Fig. 11 (a), the analysis 
made allowance for the Gaussian probabilities appropriate to the CDF determination of 
jet energies and excluded events which did not satisfy the CDF cuts as specified in the last 
paragraph. The reader should note that this is not a distribution of the mean peak mass, but 
represents rather the form of a "typical m distribution" under the circumstances specified. 
Its peak is at 137 GeV with a width (FWHM) of about 18 GeV. For Fig. 11(b) the jet 
energies input were smeared by making random changes in their magnitudes, with a standard 
deviation taken from the CDF work, intended to represent the effects of jet fragmentation 
and soft-gluon bremsstrahlung and of the detector efficiencies. This has depressed the mean 
probability distribution by a factor of about 2/3; the peak location is about 1 GeV lower 
and the FWHM increased to about 25 GeV. 

(b) ISAJET [24] 

This provides a more sophisticated simulation of tt production. A full simulation of the 

CDF detector effects has been carried out, using the appropriate jet-finding procedures and 
fragmentation codes. These detector characteristics were taken account of in the code QFL, 
as understood through the analysis of the 1988-89 CDF run [24]. The standard ISAJET 
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program does not take account of t decay to bW^, the process which dominates over t-jet 
fragmentation for the top mass values of interest to us here; however, ISAJET does have the 
option of adding one additional channel for transitions of the top quark, beyond those for 
the usual jet development, and this has been made use of here. This channel then becomes 
the dominant one, so that the later jet development is entirely that of the secondary b- 
quark from top decay and of the tertiary (c and s or -u and d) quark decays. The 500 tt 
production and decay events generated in this way for rrit = 140 GeV have been processed 
as discussed above for systems with the + Ajets) final state structure and the resulting 
mean probabihty distribution has been plotted in Fig. 12. This distribution does show a 
quite definite peak at 130 GeV, a shift of 10 GeV below the input mass, with a FWHM 
of about 25 GeV if the secondary peak at about 108 GeV is ignored. This relatively large 
shift is believed to be due the gluon radiation which is taken into account in ISAJET, but 
this is still under investigation. The secondary peak is believed to arise from the "wrong" 
combinations of jets but this is not yet certain; however, it appears unlikely to interfere with 
the use of the upper peak to determine the top quark mass, since it may be masked by the 
low-mass contributions to P{m \ event) from background events. We note also that the mass 
distributions predicted for the ISAJET model spread to much higher top-mass values than 
do the simpler "toy model" calculations. 

These model calculations demonstrate that when sets of input data obtained from 
"events" generated by an algorithm representing any of the processes 



P + P^{'^_^\}b + b+C^^''''} + X, (4.10) 



or their charge conjugates, which have proceeded through the hadronic it production pro- 
cess (2.4), are considered, the method of analysis we have proposed [23] does lead to m 
distributions which peak strongly at a mass close to the top mass used in generating these 
events, at least on the average. This holds true even when we do not know the identity of 
the quark or gluon which has generated each of the observed jets. However, we must next 
enquire whether a systematic peaking of the m distributions computed from real data at 
some particular mass value, say m*, necessarily implies the existence of a top quark with 
a mass close to m*. Are there other processes leading to the final states (/ -|- ijets), which 
we might describe as "background" , for which our analysis leads to such peaks even though 
there may be no top quark at all or perhaps just no top quark in the mass range explored? 
We shall see that the answer is "yes" because the distribution P(m) from the background 
must vanish at the threshold (Myi/ + m^) and approach zero again as m becomes sufficiently 
large since the parton parameters a and a cannot exceed unity. Precisely where its maxi- 
mum occurs depends on the nature of the background and on what cuts are made to reduce 
its contribution to P(m) in the mass region accessible for the top quark. 

4.4 Simulation of non-top background and its analysis as if top-antitop 
It is apparent that final states of the form 

p + p^l + Ajets + i^i + X (4.11) 



can readily arise without top production, through single W production and decay, as follows: 
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p + p —>■ W + {4: hadronic jets = qqqq or qqgg or gggg), W ^ I + ui. (4.12) 

Berends,et a/. [26] have provided full tree- level calculations of the cross sections for all of the 
Standard Model processes which involve one W-boson and n < 6 quarks, antiquarks and 
gluons. The case n — 6 includes all of these processes (4.11) and therefore provides a natural 
and plausible "I + 4jets" background to the "/ + Ajets" states (4.1) which result from the 
sequence of tt production (2.3), their decays (2.6) and finally leptonic decay for one W and 
decay to three quarks for the other, as outlined in (4.1). In the latter sequence, all four 
jets are quark or antiquark jets, one being a b jet and another a b jet. The calculation by 
Berends,ei al. for the processes (4.12) at the Tevatron energy predicts, after application 
of the cuts specified in Sec. 4.3 above, that the set of (one quark, one antiquark and two 
gluon)-jets occur more often (56%) than the set of (two quarks and two antiquarks)-jets 
(42%), while the set of four gluon jets is quite rare (2%). It predicts also that, with these 
cuts, b and b jets should occur among these quark jets only rarely (3%). Berends, et al. 
also calculated the rate for the tt production process (2.4) in the same approximation, with 
the decay of the consequent or W~ boson leading to the final states (4.10). Applying 
the same cuts to these final states, they then concluded from their calculation that, for top 
masses between 100 and 135 GeV, the number of these events due to it production and decay 
would exceed the number of background events. 

The m distributions obtained by subjecting these QCD background events to our 
analysis procedure may be expected to differ appreciably from those for real tt production 
and decay. We have examined this question by using the VECBOS Monte Carlo program[27] 
which implements these Standard Model calculations of Berends, et al. for "W+3jets" and 
"W-|-4jets" events, to obtain a large sample of calculated events which are input to the 
ISAJET+QFL evolution and development. The sample of "W+Sjets" used corresponds to 
an IL = 112 pb~^ and that of "W+4jets" to IL = 128 pb~^. Each event was put through 
our analysis procedure as described above, assuming it to result from tt production and 
decay, and gave an m probability distribution P{m \ I + 4jet). The mean m distribution, 
normalised to JL = 4 pb~^, is shown by the shaded area on Fig. 13. It does show a peak 
at m 112 GeV, with an appreciable tail running up as high as 140 GeV. The open area 
of Fig. 13 corresponds to the fictitious event Fi shown on Fig.8, with its peak at 135 GeV; 
the two peaks are quite well separated, showing that a tt event of type Fi with top mass 
135 GeV could readily be separated from this background, given sufficient events. However, 
if the top mass were much lower than 120 GeV - the value our analysis of the data on CDF-1 
might suggest - the "l + Ajets" m distribution would overlap so strongly with the background 
m distribution as to make their separation more problematical, requiring a very large body 
of data, at the least. In this situation, the most convincing evidence on rrit would be that 
obtainable from the (/x^e"^) dilepton events, where no source of serious background is known. 

4.5 b— and/or 6— quark tagging. 

This situation could be much improved if we had some knowledge of the identities 
of the four jets in these "/ + ijets" events, since some permutations of the labels (1,2,3,4) 
would then be excluded and the spread of the P{mt \ I + Ajets) distribution due to incorrect 
identifications would be much reduced. For example, for our fictitious event F^e; if jet(0) 
had an accompanying low-energy lepton in roughly the same direction, this would suggest 



20 



that jet(0) in Table 4 is the b jet, since the two leptons e~ and /"*" would then stem naturally 
from the i decay sequence t bW^, followed by 6 ^ cl^i^i and W~ e~i?e- In this case, 
jet(0) would be assigned the label 4 and the only permutations which need be made are those 
within (1,2,3). On the other hand, if the secondary lepton were an l~ , then jet(0) would be 
the b jet and should be assigned the label 3 in the above discussion; only the permutations 
within (1,2,4) would be necessary in the sum for the net probability. In Sec. 4.1, we noted 
that the emission of a lepton from either the b or the b quark will occur in about 42% of the 
it events (2.5), so that observation of one secondary lepton should be a common occurrence. 
Sometimes this lepton may escape detection because of its relatively low energy, so that, in 
practice, its rate may not be as high as this remark would suggest, but whenever it can be 
observed and measured, it can be made use of. 

Now that microvertex detectors are being used by CDF, it will often be possible to 
identify the b and b jets from the finite visible (or inferred) path up to the decay of the 
and S°(5°) mesons to which they give rise. With a magnetic field present, measurement 
of the charge signs of the B decay products will enable and B~ to be separated. The 
existence of (B^, B^) mixing, with quite a high rate, will complicate their assignment to b or 
b jet. However, the formalism to take this mixing into account is well-known and without 
free parameters (apart from those concerning CP violation, which makes no appreciable 
contribution to the hadronic phenomena under discussion here) , so that it can be built into 
the calculation of P(m | I + 4:j) from the data. Observations with the microvertex detectors 
should both support and supplement the observations of secondary leptons, but even together 
they will not make the identification of the 6 jet or 6 jet possible in all events. Nevertheless, 
such "b-tagging", whether from secondary leptons or from determining the fiight time from 
source to decay for the b quark, will prove to be of great value both for reducing the number 
of misinterpreted background events and for identifying some of the quark and antiquark 
jets, thus reducing the number of irrelevant jet assignments included in the determination 
of P(m I / + 4j) from the data. 

5 DISCUSSION 

It is clear that we can do little more at present beyond discussing possible procedures 
for the analysis of top-antitop candidate events. There is little data and it suffers serious 
imcertainties, so that we are not able even to test our theoretical assumptions yet. Never- 
theless, it is natural to assume that the top quark exists with some mass value yet to be 
determined, and that it does decay dominantly through t bW'^, as the Standard Model 
indicates for the mass values we must consider today. We beheve that we should launch 
into the discussion of all candidate events. Even if it is not possible to prove that individual 
events are necessarily top-antitop, we may still find that there is a substantial fraction of 
the {ir^e^2jets) events that give P{m \ event) distributions which are compatible and can 
be combined, in a first step to the determination of the top mass. We have illustrated this 
here by combining the P{m \ CDF — 1) distribution with the P{m \ N^g) determined from 
the rate of observed top-antitop pair production at the Tevatron; they are compatible and 
combining them gives an improved estimate of irit. 

We have noted above that "/-|-4jets" events from the decay of a it pair are expected to 
occur twelve times as often as "/i='=e^2jets" . If the three events of the latter type are really 
due to it, there should also be, in the recent data, roughly 12 x 3 = 36 events of the former 
type. At present, there do not appear to be even as many candidate events as this remark 
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suggests, whereas the majority of such candidate events are not compatible with top mass 
as high as the analysis of CDF-1 and the (^^e^) event rate would suggest. This channel 
calls for open discussion of all "/ + 4jets" events. It could prove to be a very fruitful line 
of investigation, where the advantage of much higher statistics has to be balanced against 
the necessity for discriminating against background events. It is worth remarking here that 
candidate events are a most powerful stimulus to a theoretician's thinking; it is not good 
for the progress of particle physics that experimenters should hold the details of new events 
unnecessarily long. 

Berends et al. have provided an excellent QCD model for non-top production of 
"/ + ijets" events, whose implications must continue to be studied further. Nevertheless, 
this is only a tree-level model and the proposal and study of other, no doubt more elaborate, 
models of background will merit much more attention. However, these models will not come 
forth in the absence of a knowledge of the nature of the events which are actually occurring. 
For example, there are certainly events which have more than the minimal number of jets, 
such as "l+njets" coming from tt production as well as from its background of W-production 
and decay, where n > 5. We shall need quantitative estimates of the rate of these more 
complex final states, which compete with the simpler states having the minimal number of 
jets and which so reduce the rate of the latter. 

The observation of the rate N^e of the production of the states (//^e^) for high energy 
leptons is a direct measure of tt production. The competing source of (/i^e^) is VF"*"!^" 
pair production, as pointed out in Sec. 2. 2, which is anticipated to have a lower rate. Both 
of these processes, {it + jets) and {W^W~ + jets) are also calculable at tree level. The 
former modifies the cross section for {/I'^e^), while the latter provides background. The 
determination of from the rate of (//^e^) events also needs attention to theory beyond 
tree level. The simplest calculations indicate that the N^e rate at energy 2Ep falls rapidly 
with increasing rrit, at the Tevatron energy, and it is reasonable to expect that the corrections 
just mentioned will not change this. The conclusion is that the value of rrit obtained from Nfj_e 
should be relatively stable with respect to experimental errors in N^e, such as the inclusion 
of uncertain background events, as we noted in Sec.2.2. 

6 CONCLUSION 

No firm conclusion can be expected from this work so far, except that we should be 
active in analysing all of the candidate events as they emerge, and that the " ii^e^2jets" 
events have a special value since they suffer little background from non-top events. It will 
be possible to combine the P{mt) probabilities from different events and event types as the 
data improves. One advantage of the relatively low energy of the Tevatron is that two it 
pairs will be produced in the same event only very rarely, so that this source of confusion 
will be essentially absent. Of course, a somewhat higher energy would have the advantage 
of increasing single it production beyond the low rate possible from the present Tevatron. 
The proposed increase of its cm. energy from 1800 to 2000 GeV this Summer, with other 
upgrades, will raise the rate by a factor of two; clearly an optimum energy would lie still 
somewhat higher than this. We shall not try here to estimate what would be the optimum 
energy, by balancing a high rate of single it pairs against a higher rate of multiple it pairs, 
since it would be a rather academic exercise at the moment. However, it is a question worthy 
of serious consideration to give a well-judged answer. 

Finally, we must emphasize that no top candidate yet reported has been demonstrated 
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to represent it production and decay, and that it may still be some considerable time before 
we reach that stage of certainty. However, we understand that the present plans at Fermilab 
are to run on until the end of May 1993, by which time the integrated luminosity IL should 
reach about SOpfc"^. The upgrade to 2000 GeV will take place in Summer 1993, and Tevatron 
running should begin again in the following October and run steadily all through 1994. 
Thus, it is hoped that the net IL reached by the end of 1994 may be about 200p6^^, a very 
substantial improvement over the IL up to the end of 1992, with which this paper has been 
concerned. 
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9 FIGURE CAPTIONS 

1. The partial hfetime given by the Standard Model for t — > be'^Ug is plotted vs. top 
mass rrit. The factor 9 ensures that for nit > lOOGeV the quantity plotted is the total top 
quark lifetime. 

2. The mean number of (/x^e^) events expected for integrated luminosity 30pb~^ for pp 
collisions at total energy 1800 GeV, as calculated by Eichten[15] and Crane[16]. 

3. Shows the probability distribution for m^, given that there are 3 /i^e^ events observed 
for integrated luminosity 30pb'^. 

4. Shows the eUipse PQ comprising all configurations for the vectors t,b and [ when the 
top mass is m^. 

5. The t and i ellipses for the event CDF-1, as function of m*, specified at 10 GeV 
intervals. 

6. Projections of the t and t eUipses for the CDF-1 event onto the plane perpendicular 
to the beam direction for a number of rrit values. 

7. Shows (a) P{mt \ CDF — 1) for zero initial transverse momentum and for Gaussian 
initial transverse momentum with s.d. p = 0.025mt and O.lm^; for the dashed curve, see 
text, (b) P{mt I iV^e = 3 and CDF - 1), for p = O.lm*. 

8. Shows P{nit I F) for the fictitious events F^g and F/, described in the text. 

9. The leptons and jets in a typical "/ + Ajets" configuration (actually for the fictitious 
event Fi), following tt production and decay. A low- energy secondary muon, from b or b 
decay, is added, as indicated by a short dashed line. 

10. Three randomly chosen "I + Ajets" events from the Toy Model (Sec.4.3). 

11. The mean probability distribution obtained (a) from 1000 Toy Model events and 
(b) from the same events after "smearing" (see text above, also ref.[23]) to simulate real 
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events, random changes being made in the jet energy magnitudes to represent effects of jet 
fragmentation, soft gluon bremsstrahlung, and detector efficiencies. 

12. The mean rrit distribution from our analysis of 500 production and decay events gen- 
erated by the ISAJET+QFL programs, taking into account fully the CDF detector charac- 
teristics. 

13. The mean rrit distribution from our analysis of Standard Model "Z + 4jets" events gen- 
erated by VECBOS simulating the tree-level calculations of backgrounds events by Berends 
et a/. [26]; cuts as in Sec. 4. 2. The sample was about 60 times that appropriate for the 1988-89 
CDF run. 
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